home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / TQLRAT.z / TQLRAT
Text File  |  1996-03-14  |  3KB  |  67 lines
  1.  
  2.  
  3.  
  4. ____TTTTQQQQLLLLRRRRAAAATTTT((((3333FFFF))))                                                        ____TTTTQQQQLLLLRRRRAAAATTTT((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      TQLRAT, STQLRAT  -  EISPACK routine.  This subroutine finds the
  10.      eigenvalues of a SYMMETRIC TRIDIAGONAL matrix by the rational QL method.
  11.  
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  14.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee  ttttqqqqllllrrrraaaatttt((((nnnn,,,, dddd,,,, eeee2222,,,, iiiieeeerrrrrrrr))))
  15.           iiiinnnntttteeeeggggeeeerrrr          nnnn,,,, iiiieeeerrrrrrrr
  16.           ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn dddd((((nnnn)))),,,, eeee2222((((nnnn))))
  17.  
  18.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee ssssttttqqqqllllrrrraaaatttt((((nnnn,,,, dddd,,,, eeee2222,,,, iiiieeeerrrrrrrr))))
  19.           iiiinnnntttteeeeggggeeeerrrr          nnnn,,,, iiiieeeerrrrrrrr
  20.           rrrreeeeaaaallll             dddd((((nnnn)))),,,, eeee2222((((nnnn))))
  21.  
  22.  
  23.  
  24. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  25.      On Input
  26.  
  27.      NNNN is the order of the matrix.
  28.  
  29.      DDDD contains the diagonal elements of the input matrix.
  30.  
  31.      EEEE2222 contains the squares of the subdiagonal elements of the input matrix
  32.      in its last N-1 positions.  E2(1) is arbitrary.  On Output
  33.  
  34.      DDDD contains the eigenvalues in ascending order.  If an error exit is made,
  35.      the eigenvalues are correct and ordered for indices 1,2,...IERR-1, but
  36.      may not be the smallest eigenvalues.
  37.  
  38.      EEEE2222 has been destroyed.
  39.  
  40.      IIIIEEEERRRRRRRR is set to Zero       for normal return, J          if the J-th
  41.      eigenvalue has not been
  42.         determined after 30 iterations.  Calls PYTHAG(A,B) for sqrt(A**2 +
  43.      B**2).  Questions and comments should be directed to B. S. Garbow,
  44.      APPLIED MATHEMATICS DIVISION, ARGONNE NATIONAL LABORATORY
  45.  
  46.  
  47.  
  48.  
  49.  
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.